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1. Motivation 

Recently there has been considerable interest in dimensionally reduced Yang-Mills 
theories as a means to obtain non-perturbative information on superstring theory and 
M-theory. The possible relevance of these systems to quantum gravity appeared in [U| 
through the light cone quantization of the 11- dimensional super membrane. It was argued 
in |]l| that 10-dimensional SU{N) super Yang-Mills theory reduced to one dimension - 
i. e. matrix quantum mechanics - correctly quantizes the supermembrane in the large N 
limit. The very same system, first studied (without reference to applications to quantum 
gravity) in , has recently been interpreted as a non-perturbative attempt at M-theory 
0. Unfortunately, on the technical side, very little is known about this model: It is 
suspected |jl|, that a novel, intricate large limit is required, but only few concrete 
results are available. This motivates the study of a simpler system: The complete 
reduction of Yang-Mills theory to O + O dimensions. In addition, the reduction of the 
ten-dimensional susy gauge theory path integral to a matrix integral has been at the 
heart of an alternative proposal to directly define non-perturbative IIB string theory 0]. 
More generally we may study the complete reduction of D-dimensional SU{N) Yang- 
Mills theory. Then the path integral of the field theory simplifies to an integral over 
the group's Lie algebra, with a flat measure: a Yang-Mills integral. Denoting the gauge 
potential by X^ and their superpartners by we obtain 




<N--= I 11 I 11^)1 [[dKl^W 



(1) 



with the Euclidean "action" 

1 



The a priori allowed dimensions for the reduced supersymmetric gauge theory are 
D = 3,4, 6, 10 corresponding to A/" = 2, 4, 8, 16 real supersymmetries. For the bosonic, 
i. e. the non-supersymmetric case A/" = O, we omit the Grassmann variables and may 
study all dimensions D > 2. 

There are numerous further reasons for being interested in the integrals eq.(0): 

• The susy integrals are crucial for the computation of the Witten index of the above 
mentioned quantum mechanical gauge theories, as they contribute to the so-called bulk 
part of the index (cf also [§] for an earlier calculation) . 

• In the maximally supersymmetric case, the system describes to leading order the 
statistical distribution of a system of D-instantons (or "— 1-branes"). 

• The A/" = 16 integral appears in very recent work developing a multi-instanton 
calculus for the Af = 4, D = 4 SU{oo) conformal gauge theory 0, and again in the 
large N limit of Sp(A^) and SO(A^) Af = 4 susy gauge theory 0. 

• Finally we can regard the integrals eq. (P as a version of the Eguchi-Kawai reduced 
gauge theory. The original work ^ focused on a lattice formulation and employed 
unitary matrices, while the above integrals use the hermitean gauge connections X^. 
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This is similar to |T^ ; however, we apply neither gauge fixing nor quenching prescriptions 
to the above integrals. The interesting question is whether the models eq.(P encode 
universal information on the f uli gauge field theory as — oo. 

The integrals of eq.(0) appear to be singular due to the "valleys" of the action, 
i. e. the directions in the configuration space of the X^ where all D matrices commute. 
Recent work has however proven this intuition wrong: Yang-Mills integrals do exist in 
many interesting cases. 



2. Partition functions: Convergence properties and some exact results 



Indeed, the rigorous results of p (see also p) show that for the gauge group SU (2) the 
susy integrals converge in dimensions D = 4,6, 10. The calculations are easily repeated 



for the bosonic case [|lT|, and the convergence condition D > 5 is found. Unfortunately, 
no rigorous methods exist to date for higher rank gauge groups > 3. 



In [12],|13| we developed methods to numerically test convergence of singular 



multidimensional integrals. The idea is to perform a Metropolis random walk weighted 
by the integrand, and to merely measure the autocorrelation function of subsequent 
configurations. In this approach, a unit autocorrelation function signals the presence of 
a nonintegrable singularity. 

As an illustration we plot in Fig.l the autocorrelation function of the SU (2) bosonic 
integrals. 




Fig. 1 Autocorrelation functions versus Monte Carlo time for the SU (2) bosonic integral 
with, from the left, D = 3,4,5. 

We are clearly able to reproduce the convergence condition D > 5: For D = 5 
the configurations decorrelate well and the whole integration space is properly sampled. 
(One observes increasingly improved decorrelation for D = 6,7, . . ., not shown in Fig.l.) 
In contrast for D = 3 the system gets quickly trapped in a singular configuration: The 
Markov chain gets lost in a valley, and the integral is divergent. D = 4 shows marginal 
divergence, which agrees with the exact analytical results. 
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Applying the same method to higher rank bosonic models, and to the 
supersymmetric models, we are able to map out the convergence conditions for the 
Yang-Mills integrals. They read 



D 



= 4,6,10 


and 


N>2 \ 


D = 3 


and 


iV > 4 ' 


D = A 


and 


A^>3 > 


D>5 


and 


N>2 



for Af>0 



for 7V = 



(3) 



In particular the D = 3 susy integral is divergent (see [jTT| , , [|T^ for a more detailed 
discussion of this point). 

It would be nice to have a rigorous mathematical proof of the conditions @. Some 
understanding may be gained by considering one-loop perturbative estimates of the 
integrals eq.(|l|). One has for the supersymmetric case [14 

N, D 



n 



1,11 



X, 



X E n 

G: majcimal tree (i j): link of G 



+ . . . 



(4) 



where the x^ are the diagonal components of the matrices while it is found in 
that the bosonic integrals are approximated by 

. N, D 



1,11 



yields precisely eq.( 



(5) 



However, it 



Powercounting for large separations (a;* — x- 
should be stressed that this does not prove the convergence conditions, since one has 
to worry about configurations where some separations are small, so that the one-loop 
approximation becomes invalid. 

In the supersymmetric case the value for the integrals is believed to be known: 



Z- 



■M 
D,N 



N(N+1) N-1 
2 2 TT 2 



1 

Ar2 



D 



W D = 6, 



Af 



(6) 



16 



Y.m\N — 10, 

The A/" = 16 expression was conjectured by Green and Gutperle based on a 
calculation of the D-instanton effective action of the superstring. A derivation of the 
terms to the right of the curly bracket in eq.(|^), based on cohomological deformation 
techniques, was given in [|r^. This calculation has still an important loophole (see 
comments in |ll3|). However, the formula eq.(P) was numerically checked by Monte 
Carlo techniques in ||lT|,|[T^ up to A^ ~ 5. 
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For the bosonic case no exact value of the partition function is known except for 
SU (2) where the result reads 



It would be exciting to find the generahzation of this result to higher rank gauge 
groups - this is after all the "zero-mode" contribution to the Yang-Mills partition 
function on a D-dimensional torus. 



3. Eigenvalue densities: Asymptotics and the exact SU{2) densities 

Let us shift attention from the partition functions to the correlation functions of the 
models. The simplest correlators are SU{N) invariant one-matrix correlators: the 
moments (TrX|,) of one matrix, say the D-th: X£). These are directly related to 
the distribution of eigenvalues of the matrix: If the eigenvalues of X£) are Ai, . . . , Aat, 
the eigenvalue density is defined for all as 

1 ^ 

pW = {^T.^{X-\)). (8) 

1=1 

The non-zero moments of p(A) are then given by 



1 



In an ordinary Wigner type matrix model all moments exist. Yang-Mills integrals are 
more intricate. Analytical calculations have only been performed for SU{2), and we 
found the foUowing surprising results. In the D = 4 susy integral all moments are 
infinite, even though the integral itself exists, as argued above. In the D = 6 susy 
integral the first two moments are finite and one finds 

{TTXl)n=e = l^ (TrX|,)^=6 = | (10) 

while all higher moments diverge. For the D = 10 susy integral we have exactly twelve 
finite moments which are 

297 




{TrX%)n=io = -\ - (TrX«)z, 



=10 



4096 



It would be interesting to find a geometrical or combinatorial interpretation for these 
numbers. Which densities give rise to this convergence behavior? For SU{2) we can go 
farther and find the exact densities: 

p|Jf7(A) = ^A^^(^,l 8A^) (12) 
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where [/ is the Kummer-[/ function defined as 

1 roo 

U(a, b, z) = —— / dtt''-' (1 + tf-''-^ e-'' (13) 
r(a) Jo 

Now we see that the above finiteness properties of the moments result from a rather 
curious powerhke behavior of the densities at large values of A. We have for X ^ 00 

' A-3 D = 4 



Pr^(A) 



~ < D = Q (14) 



D = 10 



This power-hke behavior is very different from Wigner type systems where the fall-off at 
infinity is at least exponential. For the D- dimensional bosonic models the density can 
be worked out as well albeit less exphcitly than in eqs. ([T^) , and one finds the asymptotic 
behavior Pd(A) ~ A'^^^. 

Moving on to higher values of A^, we are unable to analytically calculate the 
eigenvalue densities with presently known techniques. We can, however, find numerically 
exact densities using Monte Carlo methods. In Fig.2 we illustrate this by plotting the 
A^ = 2,3,4D = 4 susy half-densities (we only plot the A > O part since the densities are 
symmetric functions). In the SU {2) case the exact expression of eq.(|l^) and the Monte 
Carlo data cannot be separated on the scale of the figure. 

Now we would like to know how the SU (2) result eq.(|14D generalizes to other values 
of A^. It is impossible to extract the asymptotics from histograms such as Fig.2, since the 
tails comprise only a small number of samples. Instead, we can go back to the Markov- 
chain technique of the last section and measure the finiteness of the moments eq.(|^) 
for various values of A^. We find that for the susy integrals the qualitative behavior of 
the SU (2) case persists: In the D = 4,6, 10 integrals only the first, respectively. O, 2, 12 
moments are finite. We thus conclude, in view of eq.(0), that the asymptotic behavior 



eq. is valid for all N. This is very different from Wigner type random matrix models, 
where as N increases, the density condenses onto a compact interval. At the sharply 
defined (at A^ = 00) edge of the interval Wigner distributions show universal behavior. 
We have argued that in susy Yang-Mills matrix models no such edge exists, indicating 
that the large N physics of these models is indeed very different. 
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Fig. 2 Eigenvalue (half-)densities for susy D = 4, from the left, N = 2,3,4 

Let us furthermore compare supersymmetric and non-supersymmetric Yang-Mills 




integrals. How does the asymptotic behavior of the density eq.(|T^ change in the absence 
of susy? We have already mentioned above that for SU (2) this behavior is powerlike as 
well. Actually we can guess a general formula by looking once more at the effective one- 
loop estimate of eq.(^. For one-matrix correlators the most dangerous configuration 
stems from pulling away one coordinate from a bulk configuration of all other D — 1 
coordinates. Powercounting leads to the guess 

pz5(A) ~ A-2^(^-2H3i5-^ where N > (15) 



D-2 

The same procedure applied to the susy estimate eq.@ reproduces eq . ([T^) . We then 
verified the validity of eq.(^) by the same Monte Carlo random walk procedure as 
above, measuring the finiteness of moments. We thus notice a marked difference to the 
susy situation: In the bosonic case all moments exist as N oo for all D > 3. In 
particular we expect the eigenvalue distribution to condense onto a compact support, 
much like for Wigner type models. 



4. Wilson loops: Preliminary Results 

A further natural set of correlation functions of Yang-Mills matrix integrals are Wilson 
loops. Due to the Eguchi-Kawai mechanism 0,[]IO[, one naively expects them to 
correspond at = oo to Wilson loops in the unreduced gauge field theory. In the 
proposal of |^ for a non-perturbative definition of the IIB superstring, they have been 
interpreted as string creation operators . 



Despite the dimensional reduction of the field theory to zero dimensions we are still 
able to define an infinite set of independent Wilson loops dependent on an arbitrary 
contour C in Z)-dimensional Euclidean space: 



W(C) 



) 



(16) 
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Due to the non-commutative nature of the connections X^ and the path-ordering P, 
this is a non-trivial functional of the contour C despite the fact that the X^ are spaciaUy 
constant. In the special case of a rectangular contour with lengths L and T in the 
{UijUd) plane this simphfies to 



(17) 



We would hke to understand how the loops W behave as a function of and as a 
functional of the shape of the contours, in particular whether planar loops satisfy an 
area-law. We would also like to see whether there are any telltale differences between 
the supersymmmetric and the bosonic loops. 

In the previous sections we have shown how a numbcr of exact results may be 
derived for the simplest gauge group SU (2). We have not bcen ablc to analytically 
calculate a Wilson loop for an Yang-Mills integral even for SU (2). On the other hand 
it is possible to obtain high precision numerical results for low values of A^. In Fig.3a 
we plot the Wilson loop for a square (L = T) in the case of susy D = 4 as a function of 
L for various values of N. 



1.0 

W(R) 




20. L 40 
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0.0 
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20. 



R 



40. 



Fig. 3 D — A Susy Wilson loops (a) Sguare of side length L with N — 2,4,8 (top to 
bottom) (b) sguare, regular 16-gon and 64-gon enclosed in circle of radius R with N — A. 



The behavior for small arca is casy to understand. Indeed, for an arbitrary planar 
loop enclosing a small area of size A, it is straightforward to show, using Stokes' theorem, 
that W{C) = 1 + 2^^2^Tr[Xi,X2]2) + 0{X^). This immediately gives 

(18) 



W{C) 



1 - 



1 



AA> O 



2N D 

and agrees to high precision with the numerical data. A very curious featurc of Fig.3a 
is that the loops tend to a constant for largc area. The existcncc of this constant can be 
dcmonstratcd analytically. It should be considered a finitc A^ artifact for the following 
reasons: (a) The constant decreases with A^, as seen in Fig.Sa. (b) The constant depends 
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in various ways on the shape of the contour. We checked that by distorting the rectangle 
to a slightly irregular quadrangle the constant drops to zero for all as the size increases. 
In FigS.b we show various regular polygons approximating a circle of radius R: As we 
increase the number of edges the constant goes to zero for large R. 

It is clear from the mentioned features of the Wilson loops that there is an area 
law neither for very small nor very large areas. An intermediate region in which an area 
law holds might still be present. We checked, by going to rather larger A^, that this is 
not the case for the bosonic models. There has been an interesting suggestion that 
such an intermediate region may exist for the D = 4 supersymmetric model. 



5. Conclusions and outlook 



We showed how numerical Markov chain methods can be used to verify non-perturbative 
convergence conditions for Yang-Mills integrals with and without supersymmetry. The 
same methods may be applied to establish the convergence properties of correlation 
functions. Applying the technique to invariant correlators of a single matrix, we are 
able to accurately predict the asymptotic behavior of the eigenvalue density of Yang- 
Mills matrix models. The results demonstrate an unusual powerlaw behavior which, in 
the supersymmetric cases, persists for large A^. This indicates that the large A^ limit 
of these "new" matrix models might indeed be very different from the one of the "old" 
Wigner type models. We also demonstrated that Monte Carlo methods are capable to 
rather accurately compute various quantities relevant to these models such as partition 
functions, correlation functions, spectral distributions and Wilson loops. As opposed 
to Yang-Mills quantum mechanics [0],||T]],|^ we are confronted to a system which allows 
some non-perturbative analysis, at least for finite A^. 

Yang-Mills integrals are thus an ideal laboratory for exploring new large A^ 
techniques. Powerful analytical methods will have to be developed if we are to verify or, 
maybe more importantly, if we are to bring to good use the ideas presented in 
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